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Abstract 

A discrete version of the WKB method is developed and apphed to calculate 
the tunnel splittings between classically degenerate states of spin Hamiltoni- 
ans. The results for particular model problems are in complete accord with 
those previously found using instanton methods. The discrete WKB method 
is more elementary and also yields wavefunctions. 
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I. INTRODUCTION 



A large variety of problems in physics can be modeled in terms of a spin of large magni- 
tude which is then described by a Hamiltonian which is a polynomial in the spin components. 
It is natural in many of these cases to seek a semiclassical description. This is especially 
so when the system possesses a number of degenerate states in the classical limit, and one 
wishes to see how this degeneracy is lifted by quantum tunneling. Examples include the 
Lipkin-Meshov-Glick model in nuclear physics the rotational spectrum of polyatomic 

molecules and the Hamamoto-Mottleson model |^ for the AI = 4 staggering of the 
rotational spectra of certain superdeformed nuclei. 

Perhaps the simplest example of such a problem is provided by the Hamiltonian 

n = -kiSj + k2Sl (1.1) 

where ki > 0, ^2 > 0, and Sa (« = x, y, z) are components of the spin operator S obeying 
the usual commutation rules 

[Sa, Sp] = ieap^S^. (1.2) 

Further, the spin has a magnitude S, i.e., the operator S ■ S has eigenvalue S{S + 1). As 
S —>■ oo, the system is more and more classical, and in the limit it has two degenerate 
ground states corresponding to Sz/S = ±1. For large but finite S, the quantal analogs of 
these states will be admixed by a tunneling matrix element Ao/2. The two lowest energy 
eigenstates will then be split by an amount Aq. The goal of a tunneling calculation is to 
find an expression for Aq that is asymptotically valid as S" — oo. In analogy with the result 
for particles with a position like coordinate in a potential well, we expect this expression to 
take the form 

Ao = e-(^o+Bi+B2+-)^ (13) 

where 

5o > fii > ^2 > • • • (1.4) 

in the sense that Bi/Bq, B2/B1,... ^ as S* ^ 00. By looking at the commutation 
relations for the scaled operators S^/S, we can see that 1/5* plays the same role as h in 
conventional WKB theory, and we thus expect Bq = 0{S). As is often done for a particle 
in a well, it is convenient to call e~^" the Gamow factor, or the WKB exponent, and the 
remaining terms e~(-^i+^2+ - ) ^j^g prefactor. 

In many physical problems, it is only sensible or possible to calculate Aq to "exponential 
accuracy", i.e., to find only the leading term Bq in Eq. ( |1.3| ). As a general point of formalism, 
however, it is clearly desirable to have a proper asymptotic approximation Aq to Aq itself 
and not just to IuAq. In other words, we demand that Aq/Aq ^ 1 as S* — 00. For 
this one must know all Bn that are of 0{S^) or greater. The first correct calculation to 
such accuracy that we know of was done by Enz and Schilling who mapped the spin 
onto a canonically conjugate variable pair via a Villain transformation, and then applied an 
instanton method to the resulting problem. Great care must be taken in this approach 
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with operator ordering problems, leading to a rather complex calculation. A more natural 
approach is one based on Klauder's spin-coherent-state path integral A naive application 
of instanton methods for the tunnel splitting to this path integral |jl0[, however, only gives 
the leading S dependence correctly. In particular, terms of 0{S^) in InAg are incorrectly 
obtained [|TT],ll2[. A correct calculation within this formalism has only recently been done 



by Belinicher, Providencia, and Providencia |jT3|, who show that it is necessary to consider 



non-differentiable paths by a careful examination of the discrete time version of the path 
integral. Once again, the simplicity of the instanton approach is lost. 

It is therefore desirable to have a means of calculating the tunnel splitting using only 
elementary methods of analysis. Such an approach is provided by the discrete WKB method. 
(An alternative approach, used in Ref. [|14|, also meets the criterion of being elementary.) In 



the present context, the method consists of writing the Schrodinger equation as a recursion 
relation for the expansion coefficients of an energy eigenstate in the basis, and solving 
this recursion relation under the assumption that the coefficients in the recursion relation 
vary slowly with m, the Sz eigenvalue. By and large, the method is a straightforward 
generalization of the continuous WKB method [0, and it may also be viewed as equivalent 
to the use of semiclassical dynamics for Bloch electrons. 

The earliest appearance of this method in a physics context of which we are aware is by 



Schulten and Gordon |16], who use it to find semiclassical approximations for the 3j (Wigner) 
and 6j (Racah) symbols. A cogent review of the method in the context of atomic physics has 
been given by Braun who also cites its application in a variety of other physics settings, 
as well as earlier discussions in Russian mathematics texts. It has also been previously 
applied to spin tunneling problems by van Hemmen and Siito These authors have 

not, however, utilized the method to its full potential. Specifically, they calculate Aq for a 
class of Hamiltonians including 

n = -^Sl - aS,, (1.5) 

but only the leading S dependence is correctly found. One of the goals of this paper is to 
show that the prefactor can also be found correctly to order with little extra labor. (See 
Sec. V.) 

We will illustrate the discrete WKB method by applying it to Eqs. ( |1 . 1| ) and ( |1.5[ ). 
We will find the splittings A^, k = 0, 1,2, . . . for all low lying pairs of energy levels for 
Eq. ( |1 . 1| ) , and the lowest splitting Aq for Eq. ( |1.5| ). Our answers for Aq agree completely 
with those of Enz and Schilling [^, or Belinicher, Providencia, and Providencia |jl3[, and 



we believe they are obtained with much less effort. While the discrete WKB method does 
not have the nice geometrical structure of the spin-coherent-state path integral approach, it 
provides compensation in that it yields wavefunctions in the course of the calculation, which 
may be of further use; for evaluating matrix elements of perturbations, for example pO| , pT 



Other problems commonly solved by the continuous WKB method, such as escape from a 
metastable well, are also amenable to the method, but we do not discuss these here. 

The plan of the paper is as follows. In Sec. II we briefly review the discrete WKB 
method. Our discussion is heuristic and physically motivated. See Braun [0] for one that 
is more formal. In Sec. Ill we derive a formula [see Eq. ( |3.11| )j for the tunnel splitting for a 



symmetric double well problem analogous to that for the continuous case [^,^ . We apply 
the method to the Hamiltonian (|1.1|) in Sec. IV, and to ( p..5|) in Sec. V. We compare our 
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answers for the splittings with numerical results, and in the case of Eq. (|1.5|) to answers 
given by van Hemmen and and Siito [19|, and Scharf, Wreszinski, and van Hemmen |14 . 



II. THE DISCRETE WKB METHOD 

To introduce the discrete WKB method, let us consider the Hamiltonian and focus 

on an energy eigenfunction Denoting eigenfunctions |m) with Sz\m) = m\m) as 

usual, and writing Cm = it is easy to see that Schrodinger's equation takes the form 

of a three-term recursion relation connecting Cm-2, Cm, and Cm+2- The step Am = 2 in 
the index is inconvenient, so it is preferable to change the notation somewhat, and relabel 
the states which are connected to each other by an index j, with j = 1,2, . . . , N, where N 
equals or 5* + 1. The Schrodinger equation then takes the form 

WjCj + tjj^iCj-^i + tj,j-iCj^i = ECj, (2-1) 

where Wj = Hjj and tjj±i = TCjj±i, with tj^^^+i and tAr+i^A? being understood as zero. We 
may also assume that tjj+i 7^ for any j, since otherwise the matrix for 7i is block diagonal, 
and Eq. may be regarded as applying to one block. We can also choose the phases of 
the basis states so that all tjj±i are real. Hermiticity then yields tj+ij = t*j+i = tj,j+i- 

The above problem is equivalent to a tight-binding model for an electron on a one- 
dimensional lattice with on-site energies wj, and hopping matrix elements tjj±i, and indeed 
our notation is chosen to reflect this analogy. If Wj and tjj+i vary slowly enough with j (how 
slowly will be seen below) we may use the approximation of semiclassical electron dynamics. 
To do this we first note that for the uniform case where wj = w and tj.j+i = t for all j, the 
eigenstates are Bloch states with wavefunction e*'^-' and energy E{q) = w + 2t cosg. For the 
nonuniform case, we first extend j to be a continuous variable, and define w{j) and t{j) to 
be extensions of Wj and tjj±i via the relations 

wU) = w„ (2.2) 

which are required to hold whenever j is an integer. We further demand that w{j) and t{j) 
be smooth enough that if j/N is regarded as a quantity of order A^*^, then 

and likewise for t{j). This equation is the formal statement of 'slowly varying' coefficients. 
Whether or not it can be satisfied depends, of course, on the problem under consideration, 
and it is this condition which determines whether or not the problem is amenable to solution 
by the discrete WKB method. It is easy to see that it holds for Eq. ( |1 . 1|) whenever S ^ 1. 
We next define a local wavevector q{j) by the relation 

cosq{j) = {E-w{j))/2t{j). (2.5) 

Further, since the electron velocity is given by dE/dq in semiclassical dynamics, we define 
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v{j) = -2t{j) sin q{j). (2.6) 

Then by a straightforward repeat of the arguments employed in the continuous WKB method 
(see [|15[, e.g.), it is a simple matter to show that two linearly independent solutions to 
Eq. ( |1.1| ) are given by 

C,^-L=exp(±^ rq{j')dj'). (2.7) 

A general solution is obtained by taking linear combinations of these two solutions. 

A formal proof of Eq. ( |2.7| ) is given by Braun , along with connection formulas, Bohr- 



Sommerfeld quantization rules, etc. We will therefore limit ourselves to a few key comments 
and comparisons with the continuous WKB method. 

The first point is that Eq. (p.7|) represents a development of logCj in powers of A^~^, in 
which only the leading two terms are kept. To this accuracy, we need only keep the first 
two terms of w{j) and t{j) in an expansion in powers of 1/N. Alternatively, we can use 
any expression for these functions which will reproduce the first two terms correctly. In 
particular, Eqs. ( |2.2| ) and ( p.3|) need only hold to this order. In considering this point, it 
should be remembered that j/N is regarded as being of 0(1). 

The second point is that the factor l/^/v in Eq. ( p.7|) has the effect of normalizing each 
of the two solutions so that the probability flux is conserved and has a value of unity. This 
follows from noting that \Cj\'^ ~ 1/v, so that the particle spends a time inversely proportional 
to its velocity in any given coordinate segment dj. Since this is exactly what we expect for 
a classical particle, the quantum mechanical wavefunction is correctly adjusted to give a 
conserved probability flux. 

The third point is that for a given energy E, the classically accessible region of motion 
is defined by the inequalities 

U~ij)<E<U+{j), (2.8) 

where 

U^=wU)±2\tij)\. (2.9) 
This follows from Eqs. ( p.5|) and ( p.6|) because whenever the condition (|2.8|) is violated, q{j) 



and v{j) acquire imaginary parts. The expression ( p.TD [with yv{j) replaced by y 1^(^)1 as 
usual] then describes an exponentially decaying or growing wavefunction as opposed to an 
oscillatory one. It may also be noted that for uniform Wj and tj,j±i, the quantities f/^ would 
be the limits of the allowable range of Bloch state energies. More generally they may be 
regarded as local, j-dependent band edges. 

Finally, let us note the conditions for the validity of Eq. (p.7|) . The physical requirement 
is that the wavelength of the particle should not change by very much over dimensions of 
the wavelength itself, which is equivalent to demanding that 

|rfg(j)M1«sin2g(j)- (2.10) 

In addition to failing if w{j) and t{j) do not vary slowly (as defined above), it is easy to 
see that Eq. (|2.10|) is violated whenever g ?a or ±7r. In light of Eqs. (|2.8|) and (|2.9|) , this 
is hardly surprising, since these are turning points of the classical motion. Solutions of the 
form ( |2.7| ) on opposite sides of these points must be related to each other by the above 
mentioned connection formulas, as in the usual WKB method. 
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III. GENERAL FORMULA FOR TUNNEL SPLITTING 



In this section we consider the problem where the functions Wj and tj.j+i possess reflection 
symmetry about some point, and the central region is classically forbidden for some energies. 
The band-edge functions U^{j) then have the form shown in Fig. 1. We will derive a 
formula [see Eq. for the tunnel splitting between states formed from symmetrical 



states localized in the two potential wells. Our approach is completely analogous to that 
for a continuous coordinate |P21. Except for the extension to odd and even numbers of 



sites, essentially the same derivation also appears in Ref. p4[, and a formula very close to 
Eq. ( ^.111 ) also appears in van Hemmen and Wreszinski [^. Our notation is sufficiently 
different, and the argument sufficiently brief, that it is worth giving it in toto. 
It is convenient to shift the index j so that it runs over the values 

_N-l N-3 Jl^ (3.1) 

2 2 2 ^ ^ 

Note that j takes on integer values when N is odd, and half-integer values when N is even. 
With this labelling, we have 

tj,k — i-k,-j — ^kj- ("^-2) 

It is also convenient, for future use, to define jo = 1 or 1/2, depending on whether N is odd 
or even, respectively. In either case, jo labels the first site to the right of the mid-point. We 
also define J = (A^ — l)/2, and t^aid = ^lo or ti/2 -1/2 for odd or even A^. 

Let Cj be the normalized wavefunction localized in the right well, and which satisfies 
the Schrodinger equation neglecting the possibility of tunneling, with some energy E^. 

Further, let aj and sj be the exact associated antisymmetric and symmetric wavef unctions, 
respectively, and let them have energies Ei and E2, where {Ei+E2)/2 = Eq. These functions 
are given to very good accuracy by 



1 (3.3) 

Since the product CjC-j is everywhere exponentially small, these functions are properly 
normalized to unity. 

The Schrodinger equations obeyed by Cj and aj are 

{w, - Eo)Cj + tjj-iC,-i + t,j+iCj+i = 0, (3.4) 
{wj - Ei)aj + tjj^iaj^i + tjj+iaj+i = 0. (3.5) 

We multiply Eq. (|3.4|) by aj, Eq. (|3.5| ) by Cj, subtract the latter from the former, and sum 
the result from j = jo to J = (A^ — l)/2. This yields 

{E, - Eo) J2 Cja, + Si - S2 = 0, (3.6) 

j>jo 

where 
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El — ^ aj(tjj_iCj,_i + tjj+iCj+i), (3-7) 

j>jo 

^2 = X! + (3-8) 

i>io 

By shifting the index of summation in Eq. ( |3.8| ) we can see that most of the terms in the 

sums El and S2 are identical Recalhng that = 0, and that = for odd A^, we 
obtain 

_ f ait^idCo (oddA^); , . 

' lai/2t„.id(Ci/2 + C_v2) (eveniV). ^"^"^^ 



To evaluate the sum J2jO'jCj in Eq. (|3.6|), we use Eq. (p.3[). Since C_j is exponentially 
small everywhere in the right hand well, and since Cj is normalized, 

E C,a, - ^ E = 4k- (3.10) 



Substituting Eqs. ( p.3|) , ( p.9|) , and ( p.lO|) in Eq. ( p.6|) , and recalling that £"1 — E2 
2{Ei-Eq), we obtain |2| 



P P _/-2t„.MCo(Ci-C_i) (oddiV); 

^^"^^-j -2tn.d(C^i%-C!,/2) (even AT). (^-ll) 



Equation ( |3.11 ) is the sought for formula for the splitting. Note that unlike the case of 



a particle with a continuous coordinate degree of freedom, it is not now necessary for the 
antisymmetric state to have the higher energy. In fact, it is not hard to show that 



IV. TUNNEL SPLITTINGS FOR MODEL HAMILTONIAN (1.1) 

We now apply our general formalism to the Hamiltonian ( p..l[ ). Writing a general eigen- 
state of Ti as 

|^) = 5:C^|m), (4.1) 

m 

where Sz\m) = m\m), Schrodinger's equation becomes 
with 

= ih + \k2) {S{S + 1) - m^), (4.3) 

W+2 = \k2[{S{S + 1) - m(m + 1)){S{S + 1) - (m + l)(m + 2))]^/^ (4.4) 
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For convenience, we have added a constant kiS{S + 1) to the Hamiltonian, and thus to Wm- 
As noted earher, Eq. ([4.2|) is of the same form as Eq. (|2.1| ) with the minor difference that 
the index m jumps in steps of 2 rather than 1. We can take care of this point by changing 
the integral in the exponential in Eq. (|2.7|) to / q{m')dm' /2, leaving the other formulas 
unaffected. 

One elementary point should be noted at this stage. For half-integer S, Kramers' theorem 
ensures an exact double degeneracy of all energy eigenstates of the Hamiltonian ( p..!] ). This 
means that vanishes. In the spin-coherent-state path integral approach, this result comes 
about because the kinetic term in the action is a Berry phase, and different symmetry related 
paths acquire phases which give rise to destructive interference for half-integer S [p5|-|27[]. In 



the discrete WKB method, this result comes about from an almost trivial separation of the 
eigenvalue problem into two isomorphic disjoint subspaces. 

There are two steps in the calculation of via Eq. (|3.11| ). First, one must find the 
wavefunction Cm in the classically allowed and disallowed regions separately. In the present 
case, an expression for the former can be found exactly and in such a way that it actually 
holds in part of the disallowed region. This allows one to match on to the quasiclassical 
wavefunction given by Eq. ( |2.7| ) without using connection formulas. The second step consists 
of substituting the quasiclassical wavefunction into Eq. ( |3.11| ). 

To execute the first step, let us consider a wavefunction C^'-* with energy Eq'^ which 
corresponds to the A;th state localized in the well near m = —S, and begin by finding it 
near this well. The coefficients Wm and tm,m+2 are evaluated by writing m = —S + n, and 
expanding Eqs. ( |4.3| ) and ( [4.4| ) in powers of 1/5, with n regarded as a quantity of order 5°. 
This yields 

w^s+n= {ki + ^h)i2n + l)S, (4.5) 
t.s+n,-s+n+2 = iA;25((n + l)(n + 2)f. (4.6) 

Rewriting C^s+n = bn, the recursion relation (|4.2|) reads 



i2h + k2)S{n + \)hn + \k2S^n{n + 1)6„_2 + \k2S^J{n + l){n + 2)6„+2 = Eh^, (4.7) 
which is immediately recognizable as arising from the harmonic oscillator Hamiltonian 

Hose = (2A;i + k2)S (a^a + \) + ^k^S (a^ + (a^y) , (4.8) 

where a and are the usual annihilation and creation operators obeying [a, = 1. To see 
this let us denote the A;th energy state of Hose by li^k), and the number eigenstates by \n) 
with a'^a\n) = n\n). It is then apparent that b^^^ = {n\ipk) obeys the eigenvalue equation 

To find bl^\ we transform to position and momentum operators x and p obeying [x, p] = i 
via a = {x + ip)/ \f2, = {x — ip) / ^/2. In terms of these. 

Hose = {ki + Ik^) S{x' + p") + \k2S{x^ - p^). (4.9) 

It follows that the A;th state has an energy 
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4'^ = (A; + |Vo, (4.10) 

uJo = 2Sih{h + k2)f\ (4.11) 

The wavefunction of this state is given by a textbook formula: 

= {2'\ky)M'y"\-^'''^'H,{x/0. (4.12) 
where is the A;th Hermite polynomial, and 

It is also useful to write ^ = e~^, in terms of which we have the symmetrical relations 

= tanh(/€/2), 

e"" = tanh^. ^ ' ' 
The number states \n) on the other hand, have wavefunctions 

{x\n) = (22'^(n!)V)''^%-"'/2^„(x). (4.15) 

We can thus evaluate the overlap {n\il)k) = b^n'' ^ representation. The requisite inte- 
grals are most simply found by using the generating function for the Hermite polynomials, 

oo n 

g2..-.^ = Y: -Mx). (4.16) 

n=0 

The final result for ft^'^^ is 

1/2 [fc/2] 



- m ( V (-1)^' (tanhg)^+^^' (2sechg)^-^^' 

" -coslTT/;) l^^n+fccosh^j f^^ {D+jy.{k-2jy.j\ 



(4.17) 



Here D = {n — k)/2, and [x] denotes the largest integer less than or equal to x. Note that the 
cos(7ri5) factor accounts for the vanishing of 6^*^'' if n — A; is odd. The latter fact is a simple 
consequence of the parity relation if„(— x) = (— l)"iJ„(x). It can be seen as a restatement 
of the separation of Ti into two subspaces, or of the Am = 2 step in Eq. ( [4.2|) , and rigorously 
implies C^-* = for odd m + S — k for the entire range of m. 

Next, let us find C^^ in the central or classically forbidden region where S — m and S + m 
are both 0{S). In this region, we have to sufficient accuracy 

w{m) = [ki + ifca) {{S + - m2) + 0(S°), (4.18) 
tM = \k, ({S + I)' - m2) + 0(S°). (4.19) 



It is easy to see that d{w/t)/dm = 0{S ^), so the quasiclassicality condition ( p.lO|) holds if 
S ^ 1. Equation (Of) then yields for the local wavevector 
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cosg(mj 







1^ — m^^ 









'2h + k2\ , 2Eq 



V «^'2 / ^ h{2S + l)[s + l-m^ S + l + m) ' 

For the discrete WKB method to work, the energy levels which are split must lie well below 
the barrier, i.e., Eq'^ <C {ki + ^2/2)5'^. The second term in Eq. ( |4.2CI| ) is then of order S""^ 
relative to the first, and we can solve for g(m) as an expansion in powers of 1/S. Equation 
( [4.13| ) implies that 

l + 2A;i/A;2 = cosh/t, (4.21) 

so that 



g(m) = TT -\- i 



E'iU 1 1 \ ^( e!,'^ 



+ ^:—. +0 



^0 \S + I — m S + I + m 1 \ ujqS 



(4.22) 



The exponent in Eq. (|2.7| ) therefore equals (recalling the extra factor of 1/2 due to the step 
Am = 2) 

/■™ , ,,dm' mm nm e'^\ ( S + \ + m\ 
= = — + (4.23) 

Similarly, correct to leading order in S", we have 

Hm)\ = ^({S + ^f~m^) (4.24) 

We now substitute Eqs. ( ^.23) ) and ( ^.241 ) into Eq. ( ^.7] ), changing the normalization from 
unit flux to total unit probability in accord with the discussion in Sec. III. Writing the result 
only for the nonzero i.e., for m + S — k even, we have 

where Ak is the normalization constant, and we have used Eq^^ /uoq = A; + |. 

It is now possible to match Eqs. ( ^.25| ) and ( [4.17] ) directly, since both results are valid 
in the region S* ^ n ^ 1. We first find h^^'> foi S ^ n ^ k, with k = 0(1): 

^(,) ^ cos(7r(n-fc)/2) /2y/4 {secheY''+'y'{tanhe)(--'y'n^''-'y\ (4.26) 



To find in the overlap region, we put m = —S + n in Eq. ( [4. 251 ), recall that e = tanh 9, 
and in this way obtain for n ^ 1, 

= cos(7r(n - fc)/2)e-^/^(tanh^)("-^)/^n(^'=-i)/^ (4.27) 
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Comparing with Eq. ([4.26|) we obtain 

•2x1/4 ^ 



= ( -) ^(sech0)(2'^+i)/2(^g^^i^^)(5-fc)/2(25){2fc+3)/4_ ^4 28) 



The last step is to apply Eq. for the splitting. The only issue requiring caution is 



that of even versus odd A^, the number of points on which the wavefunction C^-* is nonzero. 
We have N = S + 1 for even k and N = S for odd k. Thus the parity of is that of 
S + k — 1. Both cases can be easily handled at the same time, however. From Eq. ( [4.25|) , we 



see that for \m\ <^ S, the nonvanishing C^^'s are given by {Ak/ S)e~'^"^^'^ up to an overall 
phase. Thus, the factor multiplying 2tmid in Eq. (|3.11|) equals 2sinhK l^fc/Sp in magnitude 
in both cases. Similarly, we see from Eq. (|4.4|) that tmid equals k2S'^ /A in both cases. Thus, 



Afc = fcasinhK (4.29) 

for all k and S. 

It remains to substitute Eq. ([4.28| ) into Eq. (|4.29| ). We quote the answer in terms of the 
parameters and 0: 

1 / 45* \'' 

/Cx 1/2 

Ao = 4oJo f - j sech^(tanh^)^. (4.31) 

For some purposes, it may be desirable to write these results in terms of the ratio ^ defined 
in Eq. ( |4.13| ). The conversion is easily done using the formulas 



tanh^^ = (l-e')/(l + e'), (4.32) 
sech^ = 2^/(1 + ^2)^ (433) 

(l/sinh2^) = 2e/(l-a- (4.34) 

As mentioned before, the result ( [4.31| ) for the ground state splitting agrees with that of 
Enz and Schilling |M] and Belinicher, Providencia, and Providencia A comparison with 



numerically obtained exact answers is given in Table I. As is evident, the agreement is quite 
good, and improves with increasing S. 

The formula ([4.3CI|) for the higher splitings is, to our knowledge, new. 



V. TUNNEL SPLITTING FOR MODEL HAMILTONIAN (1.5) 

In this section we apply our method to the Hamiltonian (1.5), limiting ourselves to finding 
the ground state tunnel sphtting Aq. This time the recursion relation connects contiguous 
coefficients C^-i, Cm, and C^+i, and 

w{m) = — 7m^, (5.1) 
t{m) = {-a/2){S'^ - m^y^\ (5.2) 

where we have introduced the convenient abbreviation 
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S' + 



(5.3) 



Let us begin, as in the previous section, by finding the wavefunction for the state locahzed 
in the right hand well, C^, in the classically accessible region. This time the answer depends 
on how the parameters a and 7 are taken to vary as 5" — > 00. To see this note that the lower 
band-edge function U~{m) now has a quadratic as opposed to linear form near its minima 
(see Fig. 2). If we define 

m = ,S'cos^, (5.4) 

then 

U- im) = -7^'^ cos^ e - aS' sin 6. (5.5) 

The mimima are located at 9 — 9i, tt — 9i, where 

sin ^1 = a/275'. (5.6) 

It follows that if a and 7 remain fixed as — > cxo, the minima approach the points m = ±S, 
which gives us a problem similar to that of the previous section. For novelty's sake, let us 
therefore take a and 7 to scale so that 9i remains fixed and not close to as 5" —> 00. We 
shall see that our final answer for Aq will hold even when this condition is not satisfied. 

For future use, let us define ^0 as the angle at which the classical expression for the 
energy (not U~{m)) has a minimum. We have 

sin^o = a;/27,5. (5.7) 

Further, 

sin 6*1 = (1 - (1/25) H ) sin 6*0; cos 6*1 = cos + sin^ 6*0/25 cos 6*0 H . (5.8) 

With the above conditions on ^1, it is easy to find Cm near the potential minima. Let 
us define 

mi = 5' cos ^1. (5.9) 

The conditions on a and 7 imply that S — mi = 0{S). Since the wavefunction we seek is 
essentially that for the ground state of a harmonic oscillator, its spatial extent. Am, on the 
other hand, is of order 5^/^. This can be seen by expanding the semiclassical Hamiltonian 
in powers of q and m — mi. We have 

Hsc{q, fn) = w{m) + 2t(m) cos q 

fti -75'2(1 + sin^ ^1) + 75'2 sin^ 9iq^ + 7 cot^ 9i{m - mif + • • • . (5.10) 

The leading neglected terms are, ignoring S' independent constants, of order 5'(m — mi)g^ = 
/ii, S''^q'^ = h2, and (m — mi)^ / S' = h^. It follows that the small oscillation frequency ujq 
and the ground state energy Eq (in the absence of tunneling, needless to say) are given by 

a;o = 275'cos^i, (5.11) 
Eo ^ --fS'\l + sin^ 9i) + -fS' cos 9i. (5.12) 
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The wavefunction can be written down directly as for a harmonic oscillator with a continuous 
spatial coordinate, 



Cn 



cos6'i 



1/4 



TiS' sin' 01 



exp 



cos6'i 



2S' sin' 01 



[m — rrii. 



(5.13) 



which is already normalized to unity inside the well. 

What is the range of validity of Eq. ( |5.13| )? This can be found by demanding that 
the terms ignored, hi, h2, and h^, be much smaller than ujq/2, the value of the and 
(m — mi)^ terms in Eq. ( |5.10|) . Since the kinetic and potential energies are equal for a 
harmonic oscillator, we have oc (m — rriiY / S'"^ in the region of validity. Thus, hi = 
0{m — niiY/S' = 0{h^), and /12 = 0{m — rriiY / S'"^ -C hi, /13. Demanding that ujq/2 > hi, 
we find that Eq. (|5.13|) is valid provided 



\m — mi\ <^ {S 



f^2/3 



(5.14) 



Note that this range keeps us clear of the end point m = S. 

The second step is to find Cm in the classically forbidden region, especially near the point 
m = (if 5* is integral) or m = 1/2 (if S" is half-integral). Since q{m) is imaginary, we define 



q[m) 



iK[m] 



with n{m) > 0. Using the basic definition ( p.5|) with the energy ( |5.12|) , we obtain 

sin^ 9 + sin^ 9i cos 9i 



COShfi;(m) ^ n n n n 

2 sin 9 sin 9i 2S' sin 9 sin 9i 

Writing 

K(m) = Ko{m) + /ti(m) + ■ • • , 
where kj is of order 1/S'\ we obtain 

_ 1 /y2_^2X 



+ 0{S 



l\~2 



Ki(m) 



2^^\S'^sm^9i 
1 
2 



_m — nil + ^1 
The exponent in Eq. ( |2.7D therefore equals 

0(m) = (poim) + 0i(m) 

= - [(5" + m) log(S" + m) — (S" — m) log(S" — m) — 2m log(eS" sin ^1)] 

+ - [log(mi — m) — log(mi + m)] . 



(5.15) 

(5.16) 
(5.17) 

(5.18) 
(5.19) 



(5.20) 



Here, 0o and (pi are the quantities given on the first and second lines above, and correspond, 
resepctively, to hq and ni. 

The velocity is now easy to find. To leading order in S, we get 
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|f(m)| = 275' sin ^1 sin X sinh = 7(^1 ^ ""^ )• (5-21) 
The wavefunction in the central region therefore equals 

Cm = ■ exp0oM) (5-22) 

m + mi 

where A is a constant that must be found by matching the solution with Eq. (|5.13| ). This in 
turn is done by expanding (j)o{m) in powers of m—mi, which is permissible as \m—mi \ <^ S"^^^ 
in the matching region. We find 

COS^l / \2 

(f)o{m) = 000 - ■ 2/) - ^i) ; 5.23 
25' sm Oi 

000 = 5'[logcot(ei/2) - cos^i]. (5.24) 

Putting m = nil in the pre-exponential factor in Eq. (|5.22| ), we see that the forbidden 
region Cm has exactly the structure required to match on to Eq. ( |5.13| ). Comparing the two 
expressions, we obtain 

A = 2mJ (5-25) 
\7iS' sm Oi J 



The final step is to apply the formula ( |3.11| ) for the splitting. In doing this, the following 
simplifications are useful. First, 

0o(^) = — mlogsin^i for m = 0(1). (5.26) 

Second, to leading order in S, 

Wo = 275 cos ^0- (5.27) 
Collecting together all these facts, we obtain, for both integral and half-integral S, 

S^'/' COs3/2^o/l-COS^^o\"^'/\25cos.o 



Ao = 4wo - (5.28) 

V TT / sm 6'o \ 1 + cos Oq J 

This is the main result of this section. It can be seen to be identical to Eq. (16) of Enz and 
Schilling |^, via the substitutions A = •j, a = sin 6^0 and b = 0+. It is almost, but not quite 
identical to Eq. (4.11) of Garg and Kim ||12|, which errs by a factor of 0{S^). 

Several comments about Eq. ( p. 2^ ) are in order. First, although we have only derived 
it assuming that sin^o ^ 1/5, it holds in fact for all sin6'o, provided of course that the 
WKB approximation itself is valid. This can be seen from the numerical comparison we 
give below, and also analytically, by assuming that sin^^o = 0(1/ S) from the very start. 
This suggests that Eq. ( |5.28| ) holds uniformly for all sin^o not too close to 1, but we have 
not investigated this point carefully. As a subcomment, note that in the limit 6q —>■ 0, i.e., 
vanishing a, Aq ~ a"^^ . This is of course, exactly what is expected from perturbation theory 
— it takes 25 applications of the perturbation —aS^ to connect the states with 5^ = ±5" 
||28|. Second, the WKB validity restriction means that the field a in Eq. (|1.5|) can not be 
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too large. Quantitatively, we require that the energy barrier be much greater than tuo, the 
level spacing in the well. This condition can be shown to be equivalent to 



1 -sin^o > 5"^^^. (5.29) 

Third, as long as Eq. ( |5.29| ) holds, in the limit where the field is large and the barrier is 
small, Eq. ( ^.28|) reduces to 

Ao = 4y3c<;o(5o/27r)^/2g-Bo^ (^5 30) 

where Bq = (2/3)5* cos'^ 6*0 = (16/3)Af//a;o, Af/ being the energy barrier. This form is 
the well known answer for the tunnel spliting of a massive particle moving in a quartic 
potential double well. The spin problem can be transformed into such a problem by noting 
that (i) in this limit Sx ~ S, so that the commutator [Sy, Sz] = iSx can be approximately 
reproduced by writing Sz = x, Sy = —Sp, where p and x are momentum and position 
operators with the commutator [p, x] = —i, and (ii) expanding the square root in the relation 

Sx = S"^ — {Sy + S^) , and rewriting Sy and Sz in terms of p and x as indicated above. 
We leave it to the reader to check that the frequency ujq and energy barrier AU of this 
particle problem are given by 276* cos and cos^ 6*0/4, so that the answer for Aq is 
indeed given by Eq. ( |5.30|) . 



Fourth, let us compare our result for Aq with the numerical and analytic answers obtained 
by van Hemmen and Siito ||19| , which apply when sin 6*0 = 0{1/S). This comparison is shown 



in Table II. The parameters are 7 = a = 1, and the two approximate forms for the splitting 
are given by 

Aasymp = 2^ S exp[-2S \ogi^ S)], (5.31) 



A ^0 

Aint = Y exp 



, _i / —Eq — 7m- 



rmt 

/ cosh 

J ~mt 



2 



^a[S{S + l)-m 



21 1/2 



(5.32) 



In Eq. ( |5.32| ), im^ are the classical turning points, where Eq = w{m), and Eq is the 
numerically obtained ground state energy, van Hemmen and Siito describe Aasymp as an 
"asymptotic" result, and Aint as the "exact WKB" answer. As can be seen from the columns 
containing the ratios of the splittings given by these formulas to the numerically found 
splittings, neither result asymptotes to the correct answer. [Recall that two functions f{x) 
and g{x) are said to be asymptotically equal as x — > Xq, i.e., f{x) ~ g{x), \i f{x) / g{x) — > 1 as 
X Xq] Aasymp actually gets worse with increasing S", which indicates that a multiplicative 
term such as e"*^*^ has been omitted (c is some constant), while Ai^t is off by a term of S^ 
p9| . Neither result therefore meets the goal set in Sec. I of our paper: to find all terms 



in Eq. ( |1.3| ) that are of order S^ or greater. Equation ( |5.28| ), on the other hand, does just 
that. It is also obvious that the recipe given by van Hemmen and Siito [see Eq. (5.19) of 
Ref. [|19|] cannot really be exact in general as regards the prefactor. 

Lastly, Eq. ( ^.28 ) can be compared with the answer obtained by Scharf, Wreszinski, and 



van Hemmen |T^, which is valid when sin^o = 0{1/S). Their result — see their Eqs. (1.7) 
and (1.8), and Table 1 — is also off by 0(5°). In fact the discrepancy is again (c/tt)^/^, 
which is to be expected in light of their method ||29||. 
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FIGURES 



FIG. 1. Discrete WKB potential energy diagram for tunnel splitting problem. For a particle 
with energy Eq, the region between ji and j2 is classically forbidden. 

FIG. 2. Sketch of the band-edge functions for the Hamiltonian (|1.5| ). The key point is that 
f7"(m) now has a quadratic minimum. 
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TABLES 



TABLE L Comparison between numerical and analytical [Eq. (4.31)] results for the ground 
state tunnel splitting for the Hamiltonian (^]^). Numbers in parentheses give the power of 10 
multiplying the answer. We have chosen ki = 5.0/ S, k2 = 20.0/ S. This scaling keeps fixed as 
S 



oo. 



S Aq (numerical) Ao (analytic) Error(%) 



10 


9.282(-3) 


9.749(-3) 


5.0 


11 


3.738(-3) 


3.906(-3) 


4.5 


12 


1.497(-3) 


1.558(-3) 


4.1 


13 


5.974(-4) 


6.195(-4) 


3.7 


14 


2.375(-4) 


2.455(-4) 


3.4 


15 


9.412(-5) 


9.708(-5) 


3.1 


16 


3.721(-5) 


3.830(-5) 


2.9 


17 


1.468(-5) 


1.508(-5) 


2.7 


18 


5.778(-6) 


5.927(-6) 


2.6 


19 


2.271(-6) 


2.326(-6) 


2.4 


20 


8.910(-7) 


9.115(-7) 


2.3 



TABLE IL Comparison between numerical and various analytical and semianalytical results 
for the ground state tunnel splitting of the Hamiltonian (|1.5|). The parameters are a = 7 = 1. As 
in Table I, x{n) denotes the number x x 10". The quantities A„, Aasympi and Aint are taken from 
Table I of van Hemmen and Siit6[19], and are, respectively, the numerically computed splitting, 
the result ( |5.31D , and the numerically evaluated integral expression ( p.32[ ). The quantity Aq is 
given by Eq. ( ^.281 ) of this work. 





Numerics 


Analyt 


Lc/semianalytic 


results from Ref. |[l9|| 


This 


work 


s 


An 


'-^asymp 


Aasymp/ An 


Aint 


Aint/ An 


Ao 


Ao/A 


3 


1.44(-3) 


4.12(-3) 


2.86 


2.01(-3) 


1.40 


1.48(-3) 


1.03 


4 


1.18(-5) 


6.10(-5) 


5.17 


1.64(-5) 


1.39 


1.20(-5) 


1.02 


5 


5.18(-8) 


5.12(-7) 


9.88 


7.25(-8) 


1.40 


5.24(-8) 


1.01 


6 


1.43(-10) 


2.76(-9) 


19.3 


2.00(-10) 


1.40 


1.44(-10) 


1.01 


7 


2.68(-13) 


1.03(-11) 


38.5 


3.77(-13) 


1.41 


2.70(-13) 


1.01 


8 


3.66(-16) 


2.84(-14) 


77.7 


5.15(-16) 


1.41 


3.68(-16) 


1.01 


9 


3.79(-19) 


6.00(-17) 


158 


5.36(-19) 


1.41 


3.81(-19) 


1.01 


10 


3.09(-22) 


1.00(-19) 


324 


4.37(-22) 


1.41 


3.10(-22) 


1.00 


11 


2.03(-25) 


1.35(-22) 


666 


2.87(-25) 


1.41 


2.03(-25) 


1.00 



'^The numbers have been given to greater accuracy and corrected in a few places. 
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